Operational approach to the Uhlmann holonomy 
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We suggest a physical interpretation of the Uhlmann amplitude of a density operator. Given this 
interpretation we propose an operational approach to obtain the Uhlmann condition for parallelity. 
This allows us to realize parallel transport along a sequence of density operators by an iterative 
preparation procedure. At the final step the resulting Uhlmann holonomy can be determined via 
interferometric measurements. 
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I. INTRODUCTION 

If a quantum system depends on a slowly varying exter- 
nal parameter, Berry [l[ showed that there is a geometric 
phase factor associated to the path that an eigenvector of 
the corresponding Hamiltonian traverses during the evo- 
lution. These geometric phase factors were later general- 
ized by Wilczek and Zee [2] to holonomies, i.e., unitary 
state changes associated with the motion of a degener- 
ate subspace of the parameter-dependent Hamiltonian. 
In view of the Berry phase and Wilczek-Zee holonomy, 
one may ask if a phase or a holonomy can be associated 
with families of mixed states. This was answered in the 
affirmative by Uhlmann Q by introducing "amplitudes" 
of density operators, and a condition for parallelity of 
amplitudes along a family of density operators (for other 
approaches to geometric phases and holonomies of mixed 
states and their relation to the Uhlmann approach, see 

Refs. 0, i, i, 0, s, s m m m d ) . 

As mentioned above, the Berry phases and the non- 
Abelian holonomies can be given a clear physical and op- 
erational interpretation in terms of the evolution caused 
by adiabatically evolving quantum systems. One may 
also consider the evolution as caused by a sequence of pro- 
jective measurements of observables with nondegenerate 
or degenerate eigenvalues, giving rise to a Berry phase 
or a non-Abelian holonomy, respectively. The physical 
interpretation of the Uhlmann amplitudes and their par- 
allel transport is less clear. One interpretation of the 
Uhlmann amplitude [H GJ, Q3, EE [3 is that it corre- 
sponds to the state vector of a purification on a combined 
system and ancilla. Here we suggest another interpreta- 
tion, where the amplitude corresponds to an "off-diagonal 
block" of a density operator with respect to two orthog- 
onal subspaces. In this framework, we address the ques- 
tion of how to obtain an explicitly operational approach 
to the Uhlmann holonomy, which, to the knowled ge o f the 
authors, has not been previously considered (i~7l.Tl8l. [i~9| . 
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The structure of the paper is as follows. In Sec. HT1 we 
give a brief introduction to the Uhlmann holonomy. In 
Sec. IIIII we introduce our interpretation of the Uhlmann 
amplitude and show that it is possible to determine 
the amplitude using an interferometric approach. Given 
the interpretation of the amplitude we consider an op- 
erational implementation of the parallelity condition in 
Sec. I IV A[ and in Sec. IIVBI we use the parallelity con- 
dition to establish the parallel transport. In Sec. [V] we 
present a technique to generate the states needed in the 
parallel transport procedure. We generalize the approach 
to sequences of not faithful density operators (operators 
not of full rank) and introduce a preparation procedure 
for the density operators needed in the generalized case in 
Sec.|Vll The paper ends with the conclusions in Sec. lVIII 



II. UHLMANN HOLONOMY 

Consider a sequence of density operators a%, o"2, . . . , ok 
on a Hilbert space Hi . A sequence of amplitudes of these 
states are operators W\, W2, ■ ■ ■ , Wk on Hi, such that 
a k = WkWf.. In Uhlmann's terminology [3}, a density 
operator a is faithful if its range 1Z(a) coincides with the 
whole Hilbert space, i.e., if TZ(a) = Hi (Note that what 
we refer to as a faithful operator is often referred to as an 
operator of "full rank" .) For the present we shall assume 
that all density operators are faithful, and return to the 
question of unfaithful operators in Sec. EU Using the 
polar decomposition [20] the amplitudes can be written 
Wk = y/cTkVk, where Vk is unitary. The gauge- freedom 
in the Uhlmann approach is the freedom to choose the 
unitary operators Vfc. For faithful density operators Ok 
adjacent amplitudes are parallel if and only if Wl +1 Wk > 
0. Given an initial amplitude W\ and the corresponding 
unitary operator V\, the parallelity condition uniquely 
determines the sequence of amplitudes W\, W2, ■ ■ ■ , Wk, 
and unitaries Vi, V%, . . . , Vk- The Uhlmann holonomy of 
the sequence of density operators is defined as U\jhi = 

v K vl M- 

The Uhlmann approach fits naturally within the frame- 
work of differential geometry. The amplitudes are the 
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elements of the total space of the fiber bundle, with the 
set of faithful density operators as the base manifold, and 
the set of unitary operators U(N) as the fibers. More- 
over, WW^ — a gives the projection from the total space 
down to the base manifold. Finally, given a sequence 
in the base manifold of density operators, the parallclity 
condition Wl +1 Wk > induces a unique sequence in the 
total space, leading to an element of the fiber as the re- 
sulting holonomy. 



III. INTERPRETATION OF THE UHLMANN 
AMPLITUDE 

As mentioned above, our first task is to find a phys- 
ically meaningful interpretation of the Uhlmann ampli- 
tude. We regard the density operators in the given se- 
quence as operators on a Hilbert space Hi of finite di- 
mension N. In addition, we append a single qubit with 
Hilbert space H s = Sp{|0),|l)}, with |0) and |1) or- 
thonormal, and where Sp denotes the linear span. The 
total Hilbert space we denote H = Hi ® H s ■ Note that 
H can be regarded as the state space of a single parti- 
cle in the two paths of a Mach-Zehnder interferometer, 
where Hi corresponds to the internal degrees of freedom 
(e.g., spin or polarization) of the particle and |0) and |1) 
correspond to the two paths. 

We let Q(a(°\ a^) denote the set of density operators 
p on H = Hi <8> H s such that 

(omo) = V>, (i|Hi) = V ) , (i) 

i.e., Q(<j(°\ o-W) consists of those states that have the 
prescribed "marginal states" <jq and o~\ , each found with 
probability one half. We span Q(a^°\ c' 1 ') by varying 
the "off-diagonal" operator (0|p|l). What freedom do we 
have in the choice of the operator (0|p|l)? This question 
turns out to have the following answer. 

Proposition 1. p 6 Q(a^°\ a^) if and only if there 
exists an operator V on Hi such that 

p = la^ ®\Q){Q\ + \o-V 

+iV^uV^« o|o)<i| 

+ iv^VtV^W <8> |1)(0|, (2) 

and 

VV^ < 1/. (3) 

To prove this we use the following (see Lemma 13 in 
Ref. [22j): Let A, B, and C be operators on Hi, and let 
A > 0, B > 0, and 

F = A®|0)(0| + C®|0)<l| + C t <g>|l)(0|+B®|l)<l|. (4) 



Then F is positive semidefinite if and only if 

Pk{A)CPk{b) =C, A> CB q C\ (5) 

where Pti(a) an d P"R(B) denote the projectors onto the 
ranges Tl(A) and TZ(B) of A and B, respectively. In 
Eq. the symbol B Q denotes the Moore-Penrose (MP) 
pseudo inverse [2(| of B. The reason why the MP inverse 
is used is to allow us to handle those cases when A and B 
have ranges that are proper subspaces of Hi- Note that 
when B is invertible, the MP inverse coincides with the 
ordinary inverse. 

To prove Proposition [1] we first note that if p can be 
written as in Eq. then Tr(p) = 1 and p satisfies 
Eq. If we compare Eqs. (|2|) and (j4}, we can identify 
A, B, and C, and see that they satisfy the conditions in 
Eq. ([5]). From this follows that p is positive semidefinite. 
We can thus conclude that p is a density operator and 
an element of Q(cr (0) , cr (1) ). 

Now, we wish to show the converse, i.e., if p G 
Q(cr(°\o'( 1 )) then it can be written as in Eq. ©. By 
definition it follows that we can identify A = a < -°'/2 and 
B = erW /2 in Eq. Q. Since p is positive semidefinite 
it follows that C has to satisfy the conditions in Eq. ([5]) 
and thus 

~cxW > 2Ca« V. (6) 

Define V = 2v / ^W e C\/^( T T e . From Eq. © it follows 
that V satisfies VV^ < 1. Moreover, 

l -^~)V^~) = P R{am) CP naW) = C, (7) 

where the last equality follows from Eq. ©. Thus we 
have shown that p € Q(a^°\ a^) if and only if p can be 
written as in Eq. ©. This proves Proposition [T] 

Now, consider the set of density operators Q(<r, 1//JV), 
i.e., when one of the marginal states is the maximally 
mixed state. According to Eq. |J5J| it follows that 
(0|p|l) = s/aV/{2\fN). Note that the condition in 
Eq. ([3]) allows us to choose V as an arbitrary unitary 
operator, and we thus obtain 

p = V{a,W) - |0)(0| 

® |0)(1| + ^j=W^ ® |1)(0|, 

(8) 

where W is an arbitrary Uhlmann amplitude of the 
density operator er, i.e., a — WW^ . We thus have a 
physical realization of the Uhlmann amplitude as corre- 
sponding to the off-diagonal operator (0|p|l). Note that 
Q(cr, ii/N) contains more states than those correspond- 
ing to amplitudes of a. As will be seen later, these other 
states have an important role when we consider sequences 
of density operators that are not faithful. 
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Let us note some of the differences between the above 
interpretation of the Uhlmann am plit ude and the inter- 
pretation in terms of purifications [HI, [l3], [3, EH, EH • In 
the latter, the amplitude corresponds to a pure state on 
a combination of the system and an ancilla, such that the 
density operator a is retained when the ancilla is traced 
over. Similarly as for the purification interpretation, we 
consider here an extension to a larger Hilbert space, but 
in a different manner. In the purification approach we ex- 
tend the Hilbert space as Hi (£> Hi, while in the present 
approach we extend the space as Hi®H s [23[ . Since H s is 
two-dimensional it follows that Hi&H s — Hi®Hi , i.e., 
the space we use to represent the state and its amplitude 
is isomorphic to an orthogonal sum of two copies of Hi 
(one for each path of the interferometer). With respect 
to these two subspaces the amplitude is essentially car- 
ried by the off-diagonal operator, rather than the whole 
density operator T>(a, W). In some sense the amplitude 
describes the nature of the superposition between these 
two subspaces, or equivalently, the coherence of the par- 
ticle between the two paths of the interferometer. One 
can also note from Eq. ([5]) that the total state T>((T, W) 
is mixed rather than pure in general. Note that this sub- 
space approach is closely related to the inve stig ations of 
channels and interferometry in Refs. [HI, HH Ha HH . Let 
us finally point out that the subspace approach gives a 
more compact representation than the purification ap- 
proach (if dim(7i/) > 2 and a is full rank), in the sense 
that in the former case a single qubit is added to the sys- 
tem, while for latter we have to add a whole copy of the 
original system. 

Determining the amplitude. Given a state p = 
T>(a, W) the unitary part V of the amplitude W = -foV 
can be experimentally determined by first applying onto 
p the unitary operation 



1/ 



U®\l)(l\ 



(9) 



where U is a variable unitary operator on Hi- Next, a 
Hadamard gate is applied onto H Sl followed by a mea- 
surement to determine the probability p to find the state 
This probability turns out to be 



1 

P= 2 



1 



2VN 



ReTr( v / oW t ). 



(10) 



By varying U until p is maximized we determine U = V 
uniquely, if a is faithful. Thus, V can be operationally 
defined as the unitary operator giving the largest detec- 
tion probability in this setup, indirectly determining the 
amplitude W = y/crV . 



IV. REALIZATION OF THE UHLMANN 
HOLONOMY 

A. Parallelity 

Here, we address the question of implementing the 
parallelity condition between two amplitudes. Consider 



two faithful density operators a a and a b . As mentioned 
above the corresponding amplitudes are parallel if and 
only if w£w a > 0. Let {|xfe)}fc be an arbitrary or- 
thonormal basis of Hi. We denote \xk,%) = lxfc)l x ) an d 
Px = I/ ® \ x )( x \ f° r x = 0,1. Since we use the Hilbert 
space H = Hi®H s to represent a density operator and its 
amplitude, we consider two copies of H in order to com- 
pare the amplitudes of two different density operators. 
On H <g> H we define the following unitary and Hermitian 
operator: 

Z = ^| Xfc ,0)(x/,l|®|» ; l)(Xfe I 0| 

fcZ 

+ ^lx*,l)<Xfc,0|® |xfe,0)<x/,l| 



For p, 



+P Q ® P Q +Pi <g> Pi. 
• a = V(a a ,W a ) and p b = V(a b , W b ), 
E = Tr(Zp b ® Pa ) 

= \ + ^ReTr(wX). 



(11) 



(12) 



This means that the maximal value of the real and non- 
negative quantity E is reached when W b is parallel to 
W a . 

Now we use the fact that Z is a unitary operator in 
order to test the degree of parallelity between two am- 
plitudes (see Fig.[T]). Consider an "extra" qubit e whose 
Hilbert space H e is spanned by the orthonormal basis 
{|0 e ), |le)} (not to be confused with H s and the corre- 
sponding qubit in the construction of T>(a, W)). We first 
prepare the state |0 e )(0 e | ® p b ® p a on the total Hilbert 
space H e ®H®H = H e ®Hi®H s ®Hi®H s . We apply 
a Hadamard gate on qubit e, followed by an application 
of the unitary operation 



E/z = |0 e )<0 e |®Z + |l B )<l e |®i®i, 



(13) 



i.e., an application of the unitary operation Z, condi- 
tioned on the qubit e. Finally, we apply the Hadamard 
gate on qubit e and measure the probability to find e in 
state |0 e )(0 e | [27|. This procedure results in the detection 
probability 



P 



AN 



(14) 



Thus, the probability p is maximal when W b is parallel to 
W a in the Uhlmann sense. In other words, given the state 
p a = V{a ai W a ) we prepare various states p b — V{a b , W b ) 
until we find the amplitude W b that maximizes the prob- 
ability p [lH. We have thus obtained an operational 
method to find parallel amplitudes. One may note the 
similarity between this procedure and the method intro- 
duced in [!9| to estimate the trace of products of density 
operators. 

The above approach is based on the fact that Z is 
a unitary operator and consequently corresponds to a 
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FIG. 1: (Color online) The degree of parallelity between the 
amplitudes W a and Wb of the states a a and a b , respectively, 
can be tested by applying this circuit onto the states p a = 
V(a a ,W a ) and p b = V(a b ,W b ) defined by Eq. ©. A single 
"extra" qubit is prepared in state |0 e ) and a Hadamard gate 
is applied. Conditional on the |0 e ) state of the extra qubit, 
the unitary operation Z, defined in Eq. (|f f p . is applied to 
pa® pb- After the application of a second Hadamard gate on 
the extra qubit, the degree of parallelity between W a and Wb 
can be inferred from the probability to find the extra qubit 
in state |0 e ). Parallelity is obtained when the probability is 
maximal. 



state change. As mentioned above, Z is also Hermitian 
and can thus be regarded as representing an observable. 
Thus, one may consider an alternative approach where 
given the state p a = T>(a ai W a ), we prepare states pb — 
T>(cr b , Wb) until we find the amplitude Wb that results in 
the maximal expectation value of the observable. 

It is worth to point out that the maximization which 
implements the parallelity condition relates the Uhlmann 
holonomy with various well known and closely related 
quantities. We first note that Eqs. lfT2|) and (jT4j) con- 
tain ReTj^WjfW-'a). When this is maximized over all 
amplitudes Wb of a b we find that supReTr(M / ^W a ) = 
F (cr b , a a ) , where F (a b , er a ) = Tr ^/^/oba ay /a^ is the quan- 
tum fidelity. Closely related is the transition probabil- 
ity F 2 (ab,a a ) [13], and the Bures metric d(a a ,<Jb) = 
\J1 — 2F 2 (<7 a ,tJb) [U Hl| that has been proved to be 
directly related to the Uhlmann holonomy [Ta |. 



B. Parallel transport 

The procedure to find parallel amplitudes allows us to 
obtain parallel transport. Suppose we are given a se- 
quence of operators Oj on Hi for j — 1, 2, . . . , K . We 
wish to construct a sequence pj — T>(<7j ,Wj), such that 
Wj form a parallel transported sequence of Uhlmann am- 
plitudes. Suppose moreover that pi is given (in order to 
fix the initial amplitude Wi). We can now use the fol- 
lowing iterative procedure: 

• Prepare pk- 

• Vary the preparations p = T>(<Jk+i, W) over all am- 
plitudes W of &k+i until the maximum of Tr(Zp® 
Pk) is reached. 



• Let p k+ i = p. 

After the final step K we have prepared the state px 
containing the amplitude Wk — V^Tf tAjhiVq, where U\j\a 
is the Uhlmann holonomy and V\ is the unitary part of 
the chosen initial amplitude W\. The state pk can be 
modified by applying the unitary operator 

^mod=ij®|0)(0|+V 1 ®|l)<l|, (15) 

which results in the new state 

UmodPKU^ = V{cTk,\/okUtjm), (16) 



PK = 

and hence (0|pk|1) = ^/CifC/uhl/ (^VN). Given this state 
we obtain the Uhlmann holonomy t/rjM as the unitary 
operator that yields the maximal detection probability, 
as described by Eq. (| 10[) . 

Although the iterative procedure described above is re- 
alizable in principle, it is no doubt the case that it would 
be challenging in practice, since at each step of the pro- 
cedure we must implement an optimization to find paral- 
lel amplitudes. However, the functions we optimize over 
have rather favorable properties. In the case of faithful 
density operators one can show that the function (taken 
over all amplitudes Wb) defined by Eq. (|14[) is such that 
there is no local maximum except for the global maxi- 
mum. In the case of not faithful density operators the 
global maximum is not unique, but any of them gives 
the desired result, as shown in Sec. IVII Moreover, it still 
remains the case that every local maximum is a global 
maximum. Hence, in both the faithful and unfaithful 
case we can apply local optimization methods (see, e.g., 
[33|]). The fact that local optimization techniques are ap- 
plicable is favorable for practical implementations, and 
improve the chances to find efficient procedures. How- 
ever, a more detailed analysis would be required to de- 
termine what efficiency that ultimately can be obtained. 
This question is, however, not considered here. 



V. STATE PREPARATION 

Since the parallel transport procedure involves re- 
peated preparations of states V(a,W), with arbitrary 
amplitudes W of a, we here consider preparation tech- 
niques for such states (see Fig. [5]). First, we show how to 
prepare the state p = T){a,y/o). Consider the following 
orthogonal but not normalized vectors: 



|^ = vW2|fc>|0) + |fc>|l>/V22V, 



(17) 



where Afc and |fc) are eigenvalues and correspond- 
ing orthonormal eigenvectors of a. One can check 
that J2k IV'fcKV'fcl = P- The probability distribu- 
tion (1/N, . . . ,1/N) is majorized [34| by the vector 
(Ai, . . . , Ajy)- Thus, there exists [351 ] a unitary matrix 



U such that i y,-fcl a A fc = 1/N for all j = 1,. 
(see also Refs. |13.|36|). Define the vectors 

\ m )=yfN~Yl U ^}- 



, N 



(18) 
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FIG. 2: (Color online) Preparation method to obtain the 
states p — T>(a, W) that represent density operators a and 
their amplitudes W, as defined in Eq. jB}. The output state 
p describes both the path state and the internal state of the 
particle. All states T>(a, W) can be prepared by letting a par- 
ticle in a pure internal reference state \rj} and path state |0) 
fall onto a 50-50 beam-splitter, followed by unitary operations 
acting separately in the two paths on the internal state of the 
particle. The application of the unitary operations have to be 
coordinated by a shared output of a random generator, imple- 
menting unitary operators C/j ' and Uj in respective path, 
with probability pj — 1/N. By application of a final unitary 
V in path 1 we can obtain any desired amplitude W = \J~oV . 



One can check that these vectors are normalized. Since 
U is unitary it follows that iV -1 V\ \Wj){Vj\ = P- Thus, 
p is the result if we prepare \r\j) with probability 1/N. 
One can check that (rjj\Po\Vj) = 1/2- Thus there exist 
normalized vectors \rjj), \rjj) £ Hi, such that 



V2 



(19) 



For any normalized \rf) G TLi there exist unitary op- 



erators U^ and Uj 1 ' such that Uj V '\i]) 

Uj \r]} — \rjj). The state p is prepared if we apply a 
Hadamard gate to the state |?7)|0), followed by the appli- 



(i) 



r(°)l 



\V 



and 



C°>, 



-C/, (1) ®|1)<1| 



cation of the unitary operator U, 
with probability pj — 1/N. In terms of an interferomet- 
ric approach we thus apply a pair of unitary operations 
Uj°\ U^\ one in each path of the interferometer, with 
the choice of pair based on the output of a random gener- 
ator shared between the two paths. This procedure leads 
to the output density operator p = T>(a, \/a). To obtain 
a state that corresponds to an arbitrary amplitude, i.e., 
7}(cr,y/aV) with V unitary, we just apply the unitary 
operation 1/ ® |0)(0| + V <g> |1)(1| onto p. 



ble sequences (defined below) of not faithful density op- 
erators 0. When the assumption of faithfulness is re- 
moved we have to review all the steps in the procedure. 
First, we note that Eqs. ([!]) - © are true irrespective 
of whether the involved density operators are faithful or 
not. By using the polar decompositions 



y/oje+ly/vk = Y V a k+l&ky/&k+lUk+l,k, (20) 

the Uhlmann holonomy can be reformulated as £/uhi = 
Uk,k-i ■ ■ -U2,i [2l|- If the density operators are 
not faithful then Eq. ([20]) does not determine Uk+i,k 
uniquely. However, if we require f/fe+i,fe to be a partial 
isometry [37j with initial space TZ(y/ak^/ak+i) and final 
space 1Z(^Ok+iy/o~k), then Eq. ([2"P|) uniquely determines 
Uk+i.k to be the partial isometry 



£4+1, fc — \J \fVk+lOk\f&k+i y/< 



(21) 



If the sequence of density operators is such that the fi- 
nal space of Uk+i,k matches the initial space of Uk+2.k+i, 
then we may define the Uhlmann holonomy as the partial 
isometry £/uhi = Uk,k-i ■ ■ ■ Ca,i ,3|- A sequence of den- 
sity operators that results in such matched initial and 
final spaces constitutes an "admissible ordered set" of 
density operators Q . Another way to express the condi- 
tion for an admissible sequence is 



(22) 



for k = 1,...,K- 2. 

Now we introduce some terminology and notation. We 
say that an operator W on H. is a subamplitude of a if 
WW' < g. It can be shown that W is a subamplitude if 
and only if it can be written W = y/aV, where VV^ < 1/. 
One may note that the physical interpretation we have 
constructed encompasses these subamplitudes. Given a 
density operator a and one of its subamplitudes W, we 
let V(a, W) denote the density operator in Eq. (jHJ) with 
the amplitude W replaced by the subamplitude W. One 
can see that when W is varied over all subamplitudes, 
then V{<j, W) spans all of Q(a, h/N). 

The following modified procedure results in the 
Uhlmann holonomy for an arbitrary admissible sequence 
of density operators. Let <jx,...,(jk be an admissible 
ordered sequence of density operators. Assume p\ — 
T>(ai, y/aJVi) is given, where we assume that y/aiVi is 
an amplitude (not a subamplitude). For k = 1, . . . , K—l: 

• Prepare p k = V(a k , y/akVk). 



VI. ADMISSIBLE SEQUENCES 

So far we have assumed that the density operators are 
faithful. Here we consider the generalization to admissi- 



Vary the preparation of p = T>(ak+i, ^/Uk+iV) with 
VV* < 1/ until the maximum of Tr(Zp ® pk) is 
reached. 



Let V*h 
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After the final step 



Uvu = VkV, 



(23) 



Note that we may reformulate the second step as a vari- 
ation of p over all Q(ak+i, lj/iV), and thus we vary over 
all possible subamplitudes of <t/- +1 . Note also that, by 
the very nature of the problem, the sequence of density 
operators a\, . . . , <jk is known to us. Thus, the projec- 
tors Pk(^+ iy/ok) ' that we are supposed to apply in each 
step of the preparation procedure, are also known to us. 
After the last step we "extract" the Uhlmann holonomy 
as described below. 

To outline of the proof of the modified procedure we 
first note the following fact. 

Lemma 1. Let A be an arbitrary operator on Hj. If V 
is such that it maximizes among all operators 

on Ttj that satisfies VV* < lj, then 



V = VUfA + Q, 

where Q satisfies P-jz(A)QP'ii(A'f) = Q> an( ^ w ^ ere ^7J(A) 
denotes the projector onto the orthogonal complement of 
the range TZ(A) of A. 

The following lemma is convenient for the proof of the 
modified procedure. 

Lemma 2. Let a%, a?, • • • j fjf be an admissible sequence 
of density operators, and suppose that the operator V k 
satisfies 

It follows that if V maximizes Re Tr{^/ <j k +i \Zo T k V k ) 
among all VV' < 1/, then 

%+i = Pn^rr^V = U k+1 . k V k (26) 

is uniquely determined and satisfies 

Vk+iV^ +1 ^ Ptz {vs —^- } . (27) 

Proof. We have to prove that if there exists an operator V 
that maximizes ReTr(^/<7fc + i y/o^VkV^) and is such that 
VV^ < lj, then this operator satisfies Eq. (12"6"]1 . Accord- 
ing to Lemma [1] (with A = ^/(Tk+iy^kVk ) it follows that 
we can write 



(24) 



V = \l y/0-k+ly/o^VkV k y/cfk^CJk+l V 'O 'k+l^/vkVk + Q, 

(28) 



where 



P ~ OP 



Q. 



(29) 



If we combine Eq. (|25|) with the assumption that 
the sequence is admissible, and thus P-R.(^/a^ y j a k _ 1 ) = 
^(v^y^+r): it can be shown that 

^/o'k+iVo z kVkV k t/oks/ak+x = V T fe+l cr fc\/ cr fe+ 1 - (^0) 



If Eq. (30]) is inserted into Eq. (J28J) we obtain V = 
Uk+i,kVk + Q- By the properties of the operator 
Q in Eq. {29]), together with n^^/a^y/a^Vk) = 
TZ^ak+iy/ok), it follows that V satisfies Eq. (|26|) . 

Now we have to prove that V k +i satisfies Eq. ([27]) . We 
again make use of the assumption that the sequence is 
admissible, and we find that 

Vk+iV k+1 = U k +i,kPii(^^—[)Ul +lk 

= ^+l,fc^K( v /^ v /?Fpr) C/ fe+:LA ; 



= P 



(31) 



where in the last equality we have used that the initial 
space of Uk+i,k is TZ(y/a ky /a k +i). Note that Eq. (j3Tj) 
implies that Vk+iV k+l < 1/. We have now proved that if 
there exists a maximizing operator V , then this operator 
satisfies Eqs. (|21)|) and (|2T|) . We finally have to prove 
that there actually does exist such an operator. If we 
let V = U k +\ lk y k , then the assumption of admissible 
sequences can be used to show that 

ReTr(^7TV^^^ t ) = ^^(v^HTV^fc+i.fc) 

= Tr^ ^/a k+ ~(T k ^/a k+ i, (32) 

which is the maximal value of KeTr(^/cr k -^i ^fokVkV^) un- 
der the assumption that VV' < 1/. This proves the 
lemma. □ 

Lemma [2] can be used to prove the modified procedure 
in an iterative manner. We begin with the first step 
of the procedure. Thus, we are given the state pi = 
D (<xi, •v/oTVi), where y/aiVi is an amplitude of o\, and 
hence V\V{ = -Pr.( CTi ). If we now let p = T>(a2,W), and 
vary W = ^fo^V over all subamplitudes of 02, we find 
the maximum of Tv(Zp pi) to be obtained when we 
reach the maximum of 

RbTt(W*Wi) = ReTr(V^2 y/c^ViV^). (33) 

According to Lemma [21 every maximizing V is such that 

Piz(V^V^)V = Now we let 5 = P K(V^V^)V 

and prepare the state P2 — 2? (02, \J02V2). We can repeat 
the above procedure in an iterative manner to find that 



V K = U k ,k-i • ■ ■ U 2 ,iV 1 = UtjuVx 



(34) 



Note that V\ is a partial isometry and thus may be com- 
pleted to a unitary operator V\ 38]. If we apply C/ mo d 
in Eq. (fT~5|) . but with Vi replaced by Vi, the resulting 
state is T>(o~k, ^/cif £%hi)- (Note that V\ is not unique, 

but this does not matter since V\V\ = V\V± for all such 
extensions.) Now we shall extract the Uhlmann holon- 
omy from the state T>{ok, v r o kUu^i) ■ Note that £/uhi 
is a partial isometry and that v / o r R r £/uhi in general is a 
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subamplitude of ok ■ In order to make sure that we in- 
deed extract the Uhlmann holonomy when we apply the 
procedure described in Sec.lIII| we first apply the projec- 
tion Piz(^ww ^tt) ® |0> <0 1 + ij ® |1)(1| onto the state 
f {o~k , y/&KUuh\)- By this filtering (post selection) we 
obtain a new normalized state p, for which 

(0|p|l) = AfPlZ(^^WW^[)^f°K U Vhl 

= A/"P TC ( v ^ v ^rr)\/o r fF-pR( v ^F v ^rT)^Uhi, 

where the constant A/" is a real nonnegative number, 
and where the last equality follows since 'R.(y/crK\/o~K-i) 
is the final space of J/tju- If w( 3 apply the extrac- 
tion procedure described in Sec. IIIII we find that the 
unitary operator U that gives the maximal detection 
probability is not uniquely determined. However, by 
using Lemma Q] one can show that every maximizing 
unitary operator U can be written U = ?7uiii + Q, 

where p rWv^t)° F R(^v^t) = Q - Hence ' 
Pn(^¥W^ai^T)U = Uuu is uniquely defined by this pro- 
cedure. We have thus found a modified procedure to 
obtain the Uhlmann holonomy for admissible sequences 
of density operators. 

Preparation procedures for unfaithful density operators. 
As a final note concerning the generalization to unfaith- 
ful density operators we show that the preparation proce- 
dure described in Eqs. (fT5| and (TlT)]) to obtain the states 
T>(a, W) with W an amplitude of a, can be modified to 
obtain states T>(a, W), with W an arbitrary subampli- 
tude of a. All subamplitudes W — y/aV can be reached 
via V such that VV' < lj. The set of operators V 
on Hi such that VV' < 1/, forms a convex set whose 
extreme points are the unitary operators on Tii, which 
follows from Lemma 21 in Ref. [22]. Thus, for every 
choice of V there exist probabilities fi n and unitaries 
V n% such that V — ^2 n p n V n . Hence, instead of apply- 
ing the unitary operator lj <g> |0)(0| + V ® at the 
end of the preparation procedure, we can instead apply 
ij®|0)(0|+F„(g)|l)(l| with probability This modified 
procedure results in the desired state p = T>(a, W). 

VII. CONCLUSIONS 

In conclusion, we present an interpretation of the 
Uhlmann amplitude that gives it a clear physical mean- 
ing and makes it a measurable object. In contrast to 



previous approaches where the amplitude resides in the 
total pure state of a twofold copy of the original system 
(a purification) , we suggest an alternative where the am- 
plitude is represented by the coherences of a mixed state 
on a composite system. This gives a more compact repre- 
sentation and also allows for a direct interferometric de- 
termination of the Uhlmann parallclity condition. Based 
on this, we reformulate the parallelity condition entirely 
in operational terms, which enables an implementation of 
parallel transport of amplitudes along a sequence of den- 
sity operators through an iterative procedure. At the end 
of this transport process the Uhlmann holonomy can be 
identified as a unitary mapping that gives the maximal 
detection probability in an interference experiment. 

In this paper, we consider the Uhlmann holonomy con- 
comitant to sequences of density operators, i.e., discrete 
families of density operators. However, the Uhlmann 
holonomy can also be associated to a smooth path of 
density operators Q, e.g., the time evolution of a quan- 
tum system. The parallel transport procedure discussed 
here is by its very nature iterative, but we can form suc- 
cessively refined discrete approximations of the desired 
path and obtain the Uhlmann holonomy within any non- 
zero error bound. The question is whether it is possi- 
ble to find an operational parallel transport procedure 
formulated explicitly for smoothly parameterized fami- 
lies of density operators. We hope that the framework 
we suggest in this paper may serve as a starting point 
for such an attempt. Given a family of density opera- 
tors u(s), one could consider a differential equation for 
the evolution of V(a(s), W(s)) (defined in Eq. ©) such 
that W{s) becomes the parallel transported amplitudes 
of p{s). However, it is far from clear whether such a dif- 
ferential equation could be given a reasonable physical 
and operational interpretation. 
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